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A brief review of the energy spectrum of primary cosmic rays above 10^° eV and the measure¬ 
ment techniques used to investigate the ultra high energy ones is given. This is followed by 
a discussion of the atmospheric shower profile of the highest energy event observed using the 
Fly’s Eye detector. Finally, cosmic ray simulation tools are discussed and used to investigate 
the heavy ion interpretation of the highest energy primary. The best way to contribute further 
to the understanding of this issue is by the collection of new and better data. 


1 Cosmic Ray Energy Spectrum 


The cosmic ray (CR) spectrum above 10^® eV (where the Sun’s magnetic field is no longer 
a concern) can be described by a series of power laws with the flux falling about 3 orders of 
magnitude for each decade increase in energy^ see Fig. 1. Above 10^^ eV, the flux becomes 
so low that direct measurements using sophisticated equipment on satellites or high altitude 
balloons are limited in detector area and in exposure time. Ground-based experiments with large 
apertures make such a low flux observable after a magnification effect in the upper atmosphere: 
the incident cosmic radiation interacts with atomic nuclei of the air molecules and produces 
extensive air showers which spread out over large areas. Continuously running using ingenious 
installations has raised the maximum observed primary particle energy to higher than 10^*^ evi 

While theoretical subtleties surrounding CR acceleration provide ample material for discus¬ 
sion, the debate about the origin of CRs up to lie knee (~ eV) has reached a consensus 

that they are produced in supernova explosionsB The change of the spectral index (from —2.7 
to —3.0) near the knee, presumably reflects a change in origin and the takeover of another, yet 
unclear type of source. The spectrum steepens further to —3.3 above 10i7.7 gY 
then flattens to an index of —2.7 at ~ eV (the ankle). A very widely held interpretation 

of the modulation features is that above the ankle a new population of CRs with extragalac- 
tic origin begins to dominate the more steeply falling Galactic population.^ The origin of the 
extragalactic channel is somewhat mysterious. 

CRs do not travel unhindered through intergalactic space, as there are several processes that 
can degrade the particles’ energy. In particular, the thermal photon background becomes highly 
blue shifted for ultrarelativistic protons. The reaction sequence pj —> —> Tr^p effectively de¬ 

grades the primary proton energy providing a strong constraint on themroximity of CR-sources, 
a phenomenon known as the Greisen-Zatsepin-Kuz’min (GZK) cutoffB The energy attenuation 
length of proton^ is shown in Fig. 2. A heavy nucleus undergoes photodisintegration in the 
microwave and infra-red backgrounds; as a result, iron nuclei do not survive fragmentation over 


This hypothesis is supported by AGASA 




Figure 1: Cosmic ray energy spectrum. 


comparable distances0 Ultra high energy gamma rays would travel even shorter paths due to 
pair production on radio photonscl Therefore, if the CR sources are all at cosmological distances, 
the observed spectrum must virtually end with the GZK cutoff at i? ~ 8 x 10^® eV. The spectral 
cutoff is less sharp for nearby sources (within 50 Mpc or so) tl The arrival directions of the trans- 
GZK events are distributed widely over the sky, with no plausible optical counterparts (such as 
sources in the Galactic Plane or in the Local Supercluster). Furthermore, the data are consis¬ 
tent with an isotropic ^stribution of sources in sharp constrast to the anisotropic distribution 
of light within 50 Mpc@ 

The difficulties encountered by conventional acceleration mechanisms in accelerating parti¬ 
cles to the highest observed energies have motivated suggestions that the underlying production 
mechanism could be of non-acceleration nature. Namely, charged and neutral primaries, mainly 
light mesons (pions) together with a small fraction (3%) of nucleons, might be produced at ex¬ 
tremely high energy by decay of supermassive elementary X particles (mx 1022 _ 1028 ev)0 
However, if this were the case, the observed spectrum should be dominated by gamma rays and 
neutrinos, in contrast to current observation! IlH Alternative explanations involve undiscovered 
neutral hadrons with masses above a few GeVP^ neutrinos pr^oducing nucleons and photons via 
resonant Z-production with the relic neutrino background Id or else neutrinos attaining cross 
sections in the millibarn range above the electroweak scaled If neutrinos/neutral-hadrons are 
primaries, they should point back to their sources, thereby enabling point source astronomy 
for the most energetic sources of flux at and above the GZK energy. However, the currmt 
GR-sample does not show significant angular correlation with powerful high-redshift sources Ed 
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Figure 2: Energy attenuation length of nucleons in the intergalactic medium. Notice that independently of the 
initial energy of the nucleon, the mean energy values approach to 100 EeV after a distance of approx 100 Mpc. 


2 Measurement Techniques 

There are two major techniques used in the detection of ultra high energy CRs. The first one 
is to build an array of sensors (scintillators, water Cerenkov tanks, muon detectors) spread over 
a large area. The detectors count the particle densities at any given moment, thus sampling 
the shower particles hitting the ground. By means of the ground lateral distributions of the 
shower one can deduce the direction, the energy and possibly the identity of the primary CR. 
The Akeno Giant Air Shower Array (AGASA) is the largest array so far constructed to measure 
extensive air showers. This array comprise 111 scintillation detectors each of 2.2 m^ spread 
over a grid of 100 km^ with 1 km spacingO The array detectors are connected and controlled 
through a sophisticated optical fibre network. The array also contains a number of shielded 
scintillation detectors which provide information of the muon content of the showers. The second 
technique, pioneered by the University of Utah’s Fly’s Eye detector^ consists in studying the 
longitudinal development of extensive air showers by detecting the fluorescence light produced 
by the interactions of the charged secondaries. Measurement of atmospheric fluorescence is 
possible only on clear, dark nights. Since mid-1998 a successor to the Fly’s Eye instrument, 
called Hi-Res, has started taken data at the Dugway site® In its final form it is expected to 
have a time-average aperture of 340 km^ sr at 10^® eV and 1000 km^ at 10^® eV. 

The Auger experiment has been conceived to measure the properties of cosmic rays above 
10^® eV with unprecedented statistical precision® The Observatory will consist of two sites 
(in the Northern and Southern hemispheres) designed to work in a hybrid detection mode, each 
covering an area of 3000 km^ with 1600 particle detectors overviewed by 4 fluorescence detectors. 
Surface array stations are water Gerenkov detectors (a cylindrical tank of 10 m^ top surface and 
1.2 m height, filled with filtered water and lined with a highly reflective material, the Cerenkov 
radiation is detected by 3 photomultiplier-tubes installed at the top) spaced 1.5 km from each 
other in an hexagonal grid. These stations will operate on battery-backed solar power and will 
communicate with a central station by using wireless local area network radio links. Event 
timing will be provided through global positioning system (GPS) recivers. Of course, the hybrid 
reconstruction (see Fig. 3) will be extremely valuable for energy calibration, but such “golden” 
events are expected to be less than 10% of the total event rate. 





Figure 3: Schematic representation of the operation of a hybrid air shower detector. 


It is worth mentioning that an array of ground-based detectors is ideally suited for involving 
school students and t heir t eachers in cosmic ray research. This has been proposed several times 
(for example, SCRODillii). 

3 Shower Profile of Highest Energy Fly’s Eye Event 

The Fly’s Eye observes an air shower as a nitrogen fluorescence light source whichmoves at 
the speed of light along the path of a high energy particle traversing the atmospheretS In other 
words, it directly detects the longitudinal development of the air cascade. The simulation of 
the shower evolution depends sensitively on the first few interactions, necessarily related to the 
quality ^f our “understandying” of hadronic collisions. The well known codes QGSJEtEj and 
SIBYLLC3 represent two of the best simulation tools to model hadronic interactions at the highest 
energies. The underlying idea behind sibyll is that the increase in the cross section is driven 
by the production of minijetsll The probability distribution for obtaining N jet pairs (with 
> pmin, where is a sharp threshold on the transverse momentum below which hard 
interactions are neglected) in a collision at energy ^/s is computed regarding elastic pp or pp 
scattering as a difractive shadow scattering associated with inelastic processes. The algorithms 
are tuned to reproduce the central and fragmentation regions data up to pp collider energies, and 
with no further adjustments they are extrapolated several orders of magnitude. On the other 
hand, in QGSJET the theory is formulated entirely in terms of Pomeron exchanges. The basic 
idea is to replace the soft Pomeron by a so-called “semihard Pomeron”, which is defined to be an 
ordinary soft Pomeron with the middle piece replaced by a QCD parton ladder. Thus, minijets 
will emerge as a part M the “semihard Pomeron”, which is itself the controlling mechanism for 
the whole interactions The different approaches used in both codes to model the underlying 
physics show clear differences in cross sections (see Fig. 4) and multiplicity predictions which 
increase with rising energys As can be seen in Fig. 5, for proton-induced showers, the differences 
become washed out as the shower front gets closer to the ground. 

Nevertheless, the footprints of the first hadronic collisions are still present in the longitudinal 






Figure 4: p-air cross sections of SIBYLL (dashed line), QGSJET (dots), and AIRES (solid Ipf) superimpossed on 
data obtained from collider experiments and cosmic ray experiments. E3 
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Figure 5: Comparison of muon lateral distributions at different atmospheric depths. Ep = 10 EeV.l23 
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Figure 6: Longitudinal development of 300 EeV showers from different primarypspecies together with the data of 

the highest energy event recorded by Fly’s Eye.tj 


development. In Fig. 6 we show the atmospheric shower profile of 300 EeV showers induced by a 
proton, a gold-nucleus and a dust grain|| obtained after Monte Carlo simulation with the program 
AIRES (version 2.2.1)E2 sibyll was used to generate hadronic interactions above 200 GeV. The 
simulated results are superimposed over the experimental data of the highest energy Fly’s Eye 
eventl23 It is clearly seen that the gold-nucleus longitudinal development better reproduces the 
data than a proton or a dust grain. However, this is not the case when the hadronic collisions 
are modelled with QGSJET. As can be seen in Eig. 7, in this case a medium mass nucleus likely 
fits the Fly’s Eye data. 

All in all, the primary chemical composition remains hidden by the hadronic interaction 
model. In light of this, in the next section we examine in more detail the speculative case of 
super heavy nuclei. 


4 Superheavy Nuclei 


It has been generally thought that ^®Ee is a significant end product of stellar evolution and 
higher mass nuclei are rare in the cosmic radiation. Strictly speaking, the atomic abundances of 
middle-weight (60 < A < 100) and heavy-weight {A > 100) elements are approximately 3 and 
5 orders of magnitude lower, respectively, than that of the iron groupc3 The synthesis of the 
stable super-heavy nuclides is classically ascribed to three different stellar mechanisms referred 
to as the s-, r-, and p-processes. The s-process results from the production of neutrons and their 
capture by pre-existing seed nuclei on time scales longer than most /3-decay lifetimes. There is 
observational evidence that such a kind of process is presently at work in a variety of chemically 
peculiar Red GiantsU and in special objects like EG Sagittaelil or SN1987A@ The abundance 
of well developed nuclides peaks at mass numbers A = 138 and A = 208. The neutron-rich (or 
r-nuclides) are synthesized when seed nuclei are subjected to a very intense neutron flux so that 
/3-decays near the line of stability are far too slow to compete with the neutron capture. It has 
long been thought that appropriate r-process conditions could be found in the hot (T > 10^° AT) 
and dense {p ~ 10^*^ — 10^^ g/cm^) neutron-rich (neutronized) material located behind the 

*’A very massive particle containing 10^ nucleons0 








Figure 7: Atmospheric cascade development of 300 EeV calcium-nucleus together with Fly’s Eye data. 
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outgoing shock in a type II supernova eventES Its abundance distribution peaks at ^ = 130 
and A = 195. The neutron-deficient (or p-nuclides) are 100-1000 times less abundant than the 
corresponding more neutron rich isobars, while their distribution roughly parallels the s- and 
r-nuclides abundance curve. It is quite clear that these nuclides cannot be made by neutron 
capture processes. It is generally believed that they are produced from existing seed nuclei of 
the s- or r-type by addition of protons (radiative proton captures), or by removal of neutrons 
(neutron photodisintegration). The explosion of the H-rich envelopes of type II supernovae has 
long been held responsible for the synthesis of these nuclidesa 

Putting all this together, starbursts appear as the natural sources able to produce relativistic 
super-heavy nuclei. These astrophysical environments are supposed to comprise a considerable 
population of O and Red Giant stars^ and we believe the supernovae ratS is as high as 0.2-0.3 
yr“^. It was recently put forward that within this type of galaxies, iron nuclei can be accelerated 
to extremely high energies if a two step process is invoked^ In a first stage, ions are diffusively 
accelerated up to a few PeV at single supernova shock waves in the nuclear region of the galaxyB 
Since the cosmic ray outflow is convection dominated, the typical residence time of the nuclei 
in the starburst results in t ~ 1 x 10^^ s. Thus, the total path traveled is substantially shorter 
than the mean free path (which scales as of a super-heavy nucleus. Those which are able 

to escape from the central region without suffering catastrophic interactions could be eventually 
re-accelerated to superhigh energies at the terminal shocks of galactic superwinds generated by 
the starburst. The mechanism efficiently improves as the charge number Z of the particle is 
increased. For this second step in the acceleration process, the photon field energy density drops 
to values of the order of the cosmic background radiation (we are now far from the starburst 
region). 

The dominant mechanism for energy losses in the bath of the universal cosmic radiation is 
the photodisintegration process. The disintegration rate R (in the system of reference where the 
microwave background radiation is at 2.73K) of an extremely high energy nucleus with Lorentz 
factor P, propagating through an isotropic soft photon background n(e) reads^ 


R = 


1 


Jo Jo 


( 1 ) 


2r2 Jq Jo 

where a stands for the total photon absortion cross section. Primed quantities refer to the 







Table 1: Giant dipole resonance parameters 
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Figure 8: Energy attenuation length of gold (solid) and iron (dashed) nuclei in the intergalactic medinm. 


rest frame of the nucleus. The total photon absortion cross section is characterized by a broad 
maximum, designated as the giant resonance, located at an energy of 12-20 MeV depending on 
the nucleus under consideration. For the medium and heavy nuclei, A > 50, the cross section 
can be well represented by a single, or in the case of the deformed nuclei, by the superposition 
of two Lorentzian curves of the form 


cj(e ) = fJo 


/2 ■p2 
0 


e'^ F 


(,2_,,2)2 + ,/2r2 


( 2 ) 


In order to make some estimates, hereaf ter we refer our calculations to a gold nucleus (the 
resonance parameters are listed in table I)cll Figure 8 shows the energy attenuation length of a 
gold nucleus due to interactions with the infrared and microwave backgroundsEl For comparison, 
we also show the energy attenuation length of an iron nucleus. We can conclude that despite the 
fact the abundances of superheavy nuclei are around 4 orders of magnitude lower respective to 
the iron group, the volume for the potential sources of extremely high energy superheavy nuclei 
(E ^ 300 EeV) increases substantially. 


5 Conclusions 

The energy spectrum of primary cosmic rays above 10^® eV has several compelling features. 
Perhaps the most impressive is the apparent occurrence of events above the GZK cutoff. This 
issue has been studied using simulation techniques and leads to a plausible interpretation of the 
highest energy cosmic ray event as due to the atmospheric interaction of an energetic heavy ion. 
The only reliable way to contribute further to the understanding of this issue is by the collection 
of new and better data in future experiments. 
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